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We study the thermodynamics near the generic (density-driven) superfluid-Mott-insulator transi- 
tion in the three-dimensional Bose-Hubbard model using the nonperturbative renormalization-group 
approach. At low energy the physics is controlled by the Gaussian fixed point and becomes universal. 
Thermodynamic quantities can then be expressed in terms of the universal scaling functions of the 
dilute Bose gas universality class while the microscopic physics enters only via two nonuniversal pa- 
rameters, namely the effective mass m* and the "scattering length" a* of the elementary excitations 
at the quantum critical point between the superfluid and Mott-insulating phase. A notable excep- 
tion is the condensate density in the superfluid phase which is proportional to the quasi-particle 
weight Zqp of the elementary excitations. The universal regime is defined by m*a*^T ^ 1 and 
m*a*^|^/i| ^ 1, or equivalently \n — nc\a''^ ^ 1, where dfi — fi — fic is the chemical potential 
shift from the quantum critical point {/j. = /j.c,T = 0) and n — fic the doping with respect to the 
commensurate density fic of the T = Mott insulator. We compute Zqp, m* and a* and find that 
they vary strongly with both the ratio t/U between hopping amplitude and on-site repulsion and 
the value of the (commensurate) density fic. Finally, we discuss the experimental observation of 
universality and the measurement of Zqp, m* and a* in a cold atomic gas in an optical lattice. 
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I. INTRODUCTION 

The low-temperature thermodynamics of a dilute ul- 
tracold Bose gas is well understood both theoretically 
and experimentally. The equation of state, e.g. the 
pressure P(/i, T) vs chemical potential and temperature, 
turns out to be "universal" to the extent that it depends 
only on a small number of parameters such as the mass 
m of the bosons and their 5- wave scattering length a, 
and is otherwise independent of other microscopic char- 
acteristics such as details of the atom-atom interaction 
potential. From a theoretical point of view, the thermo- 
dynamics of a dilute Bose gas is usually derived within 
a low-density expansion using mo?ii (or na^, with n the 
mean boson density) as the expansion parameter. 

Strong correlations in an ultracold Bose gas can be 
achieved by loading the gas into an optical lattice. It 
is then possible to induce a quantum phase transition 
between a superfluid ground state and a Mott insulat- 
ing phase by varying the strength of the lattice poten- 
tial [ ]. The main features of the Mott transition can 
be understood in the framework of the Bose-Hubbard 
model, which describes bosons moving in a lattice with 
an on-site repulsive interaction [ ]. 

In the vicinity of the Mott transition, there is no small 
parameter (such as density or interaction strength) that 
would allow us to derive the equation of state pertur- 
batively. Nevertheless, near the generic (density-driven) 
Mott transition, the thermodynamics of a Bose gas turns 
out to be similar to that of a dilute Bose gas up to some 
effective parameters. The origin of this similarity can be 
understood as follows. By varying the chemical potential 
from negative to positive values in a dilute Bose gas, one 
induces a (zero-temperature) quantum phase transition 
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FIG. 1. (Color online) Zero-temperature phase diagram of 
the Bose-Hubbard model on a cubic lattice showing the Mott 
insulators (MI) with density n = (vacuum) and n = 1, 
as well as the surrounding superfluid phase (SF). Point C at 
the tip of the Mott lobe shows the multicritical point where 
the transition occurs at fixed density n — 1. Away from this 
point, the transition is driven by a density change. The finite- 
temperature pressure P{iic^T) at point B is shown in Fig. 5. 
The zero-temperature pressure P{/j,0), condensate density 
no(/i, 0), superfluid stiffness ps(/i, 0), sound mode velocity 
c(/i, 0) and superfluid transition temperature Tc{/j.) along the 
dotted line AB are shown in Figs. 6-10. 



between a state with no particles (vacuum) and a super- 
fluid state with a finite density. This identifies the point 
/i = T = as a quantum critical point (QCP). Above the 
upper critical dimension (i+ = 2, the boson-boson inter- 
action is irrelevant (in the renormalization-group sense) 
and the critical behavior at the transition is mean-field 
like with a correlation-length exponent u = 1/2 and a 
dynamical critical exponent z = 2. Elementary excita- 
tions at the QCP are free bosons of mass m and their 
mutual interaction is determined by the s-wave scat- 
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tering length a in the low-energy limit. The depen- 
dence of the equation of state of a dilute Bose gas on 
m and a only is a direct consequence of the proximity 
of the QCP between the superfluid phase and the vac- 
uum: the thermodynamics is controlled by the QCP. It 
follows that thermodynamic quantities can be expressed 
in terms of universal scaling functions oi ji/T and an ef- 
fective temperature-dependent dimensionless interaction 
constant g{T) = 87rV2ma^T [ ]. This universal descrip- 
tion holds in the critical regime of the QCP defined by 
ma^|/i| and mo?T <C 1. 

The vacuum-superfiuid transition of a dilute Bose gas 
and the generic Mott transition of a Bose gas in an optical 
lattice belong to the same universality class. Both tran- 
sitions are governed by the same (Gaussian) fixed point. 
Elementary excitations at the QCP between the Mott in- 
sulator and the superfiuid phase are quasi-particles with 
effective mass m* and their mutual interaction is de- 
scribed by an effective "scattering length" a*. Near the 
QCP, thermodynamic quantities can be expressed with 
the universal scaling functions of the dilute Bose gas 
universality class and the nonuniversal parameters m* 
and a*. This conclusion is correct everywhere near the 
superfiuid-Mott-insulator transition except in the close 
vicinity of the multicritical points where the transition 
takes place at fixed (commensurate) density (Fig. 1). 

In this paper, we study the thermodynamics of the 
three-dimensional Bose-Hubbard model using a nonper- 
turbative renormalization-group (NPRG) approach [ - ]. 
In Sec. II, we derive scaling forms for various thermody- 
namic quantities (pressure, density, compressibility, con- 
densate density, superfiuid stiffness and superfiuid tran- 
sition temperature) of a dilute Bose gas and discuss the 
scaling functions in some limiting cases. The nonper- 
turbative renormalization-group approach to the Bose- 
Hubbard model is briefiy reviewed in Sec. III. In Sec. IV, 
we show that the NPRG approach enables to straightfor- 
wardly identify the elementary excitations at the QCP 
governing the generic Mott transition, and compute their 
effective mass m* and effective scattering length a* as 
well as their spectral weight Zqp. Zqp, m* and a* are 
calculated as a function of t/U . We then present various 
thermodynamic quantities obtained from a numerical so- 
lution of the NPRG equations and show that near the 
Mott transition they satisfy the scaling behavior charac- 
teristic of the dilute Bose gas universality class, except 
for the condensate density which is proportional to the 
quasi-particle weight. The experimental implications of 
our results are discussed in the Conclusion. 



II. THE DILUTE BOSE GAS UNIVERSALITY 
CLASS 

In this section we discuss in detail the dilute Bose gas 
universality class. We derive scaling forms for various 
thermodynamic quantities and compute the correspond- 
ing universal scaling functions in some limits. 



A. Universal scaling functions 

Let us consider a three-dimensional Bose gas described 
by the (Euclidean) action 



iFr{tP*[dr-y.+ 
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(1) 

where ?/^(r, r) is a complex field and r G [0, /3] an imag- 
inary time with j3 = 1/T the inverse temperature, ji 
denotes the chemical potential. The interaction is as- 
sumed to be local in space and the model is regularized 
by a ultraviolet momentum cutoff A. d = 3 and we set 
h = kB = ^ throughout the paper. 

The nature of the fi = T = QCP between the vac- 
uum and the superfiuid state can be understood from a 
RG analysis. Since the ja = ground state is the vacuum, 
there is no renormalization of the single-particle propa- 
gator and the correlation- length exponent u = 1/2, the 
anomalous dimension 77 = while the dynamical critical 
exponent z = 2. The dimensionless interaction constant 
g = 2mgK satisfies the (exact) RG equation 



ds 
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(with ^(1) = ^), when fiuctuation modes with momenta 
between A and A/s are integrated out (with a proper 
rescaling of fields, momenta and frequencies in order to 
restore the original value of the cutoff A) [ , ] . From 
Eq. (2), we obtain 
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for s > 1, 



where 



mg 



Air - 



-mgA 



(3) 
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is the s-wave scattering length which can be calculated 
from the action (1) by solving the two-body problem. 
g{s) is thus irrelevant (it vanishes for 5 ^ 00) and the 
only fixed point of Eq. (2) is ^ = in agreement with the 
fact that the upper critical dimension for the vacuum- 
superfiuid transition is (i+ = 2. 

There are two relevant perturbations about the Gaus- 
sian fixed point ja = T = g = 0: the chemical poten- 
tial jii and the temperature T, with scaling dimensions 
[fj] = Xjv and [T] = In a RG transformation, they 
transform as //(s) = s^l^ [i and T(^s) = s^T. In the criti- 
cal regime near the QCP, the pressure satisfies [8] 



(5) 



By choosing s ~ T~^l^ or s ~ and setting z = 

= 2 (with (i = 3), we can write the pressure in the 
scaling form 



(6) 
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or 



(7) 



The overall factor m^/^ comes from dimensional consider- 
ations while the factor l/(27r)^/^ is introduced for conve- 
nience [ ]. The energy-dependent effective (dimension- 
less) interaction constant ^(e) = g{s = A/^/2m\e\) is 
defined by 



(8) 



and is entirely determined by the mass m of the bosons 
and the scattering length a. T and Q and universal scal- 
ing functions characteristic of the three-dimensional di- 
lute Bose gas universality class. Equations (6) and (7) 
are valid in the critical regime near the QCP defined by 
ma^\li\ <C 1 and mo?T <C 1. Note that the interaction 
constant ^ is a dangerously irrelevant variable (in the RG 
sense) and therefore cannot be neglected: T and Q are 
singular functions of g{T) or ^(/i). Higher-order interac- 
tions, such as three-body interactions, are not considered 
here since they are irrelevant and give rise to subleading 
contributions to the pressure. 

Equations (6) and (7) imply scaling forms for other 
thermodynamic quantities. For example, the particle 
density and compressibility read 



(f,^(T)), 



(9) 



respectively, where we use the notation J^^'^'^\x^y) = 
didiT{x,y). 

For positive chemical potential, there is a superfluid 
transition at a temperature Tc. This transition corre- 
sponds to a singularity of the scaling function J^(x, y) 
when X = Xc{y). It follows that 



(10) 



with H a universal scaling function. Equation (10) im- 
plies that mo?Tc is a universal function of mo? ji. 
In the superfluid phase, using [8] 



no(/i,T) 
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with s ~ ^, one finds that the condensate density 
satisfies the scaling form 



(12) 



with X a universal function and g{ii) defined by (8). The 
superfluid density (or superfluid stiffness ps = ris/m) sat- 
isfies a similar scaling form. 



(13) 



Galilean invariance implies that the T = superfluid 
density ns(/i, 0) is equal to the fluid density n(/i, 0) and is 
therefore determined by the scaling function [Eq. (9)]. 
The sound mode velocity can be expressed in terms of the 
compressibility and the superfluid stiffness ps = Us/m 
(see, e.g., Ref. [9]), 



c(/i,T) 



/^(/i,T) ' 



(14) 



At zero-temperature, since Ug = n, the Bogoliubov sound 
mode velocity c(/i, 0) is equal to the macroscopic sound 
velocity [10]. 



B. Limiting cases 

For a three-dimensional Bose gas, the scaling func- 
tions can be obtained from perturbation theory (see Ap- 
pendix A). In this section, we discuss various limiting 
cases. 



1. Dilute classical gas 

When the chemical potential is large and negative, p < 
and \p\ ^ T, the system behaves as a dilute classical 
gas and the pressure takes the form 

P{^^,T)=[^f"T^/\-M/^, (15) 



which leads to 

T{x^ y) = for X < and \x\ 1, 

g{x,y) = x^l^e^l"" for x < and \x\ <C 1. 



(16) 



2. T = superfluid phase 

At low temperatures and positive chemical potential, 
the scaling functions can be obtained from the Bogoli- 
ubov theory [11, 1?]. At zero temperature. 
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(17) 
(18) 
(19) 
(20) 



and ps(/i,0) = ns(/i,0)/m = n(/i,0)/m. The first term 
in these equations is usually referred to as the mean- 
field result and the second-one as the Lee- Huang- Yang 
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correction [13, 14]. Equations (17-20) can be cast in the 
form (7,12,13) with 
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From Eqs. (18,19), we deduce the expression of the sound 
mode velocity (14), 



^(/i,0) 



(22) 



III. LATTICE NPRG 

In this section, we briefly review the NPRG approach 
to the Bose-Hubbard model defined on a cubic lattice [4, 
5]. The model is defined by the action 



-t (V^r>r'+C.C.)|, (27) 

where ^l^r{^) is a complex field, {r} denotes the sites of 
the lattice and (r, r') nearest-neighbor sites. U is the on- 
site repulsion and t the hopping amplitude. We take the 
lattice spacing as the unit length throughout the paper. 



to leading order in mo?ii. 



A. Scale-dependent effective action 



3. Quantum critical regime /i = 

At vanishing chemical potential, the condensate den- 
sity vanishes and the Bogoliubov theory reproduces the 
free boson result 



P(0,T) = -^^ln(l-e-''^-) 



<,V2)(^)''r-, 



(23) 



where C,{x) is the Riemann zeta function and C(5/2) 
1.3415. Equation (23) implies 

^(0,y)=C(5/2) for y^O. (24) 



4. Superfluid transition 

The Bogoliubov theory correctly describes the ground 
state and its elementary excitations but fails near the su- 
perfluid transition temperature Tc. In particular, it pre- 
dicts a first-order phase transition [ ]. The transition 
temperature can nevertheless be determined from a per- 
turbative approach in the normal phase, by considering 
the self-consistent one-loop correction to the self-energy 
(self-consistent Hartree-Fock approximation) , 



M = C(3/2) —— 

m \ ZTT 



3/2 



which leads to 



(25) 



(26) 



The NPRG is implemented by considering a family of 
models with action Sk = S -\- ASk indexed by a momen- 
tum scale k varying between a microscopic scale A down 
to [ , ]. ASk is a "regulator" term defined by 



Jo 



(28) 



where '^(q) is the Fourier transform of t/j^ and the sum 
over q runs over the first Brillouin zone ] — 7r,7r]^ of the 
reciprocal lattice. The cutoff function i?/c(q) modifies the 
bare dispersion = — 2t Xl^^i cos qi of the bosons. In the 
lattice scheme, ^A(q) is chosen such that the effective 
(bare) dispersion + Ra{^) vanishes [ ]. The action 
Sa = S -\- AS A then corresponds to the local limit of 
decoupled sites (vanishing hopping amplitude). 
In practice, we choose the cutoff function 

Rk{q) = -ZA,fc6fcSgn(tq)(l - ^q)6(l - ^q), (29) 



with A = v2(i, Ck = tP, = (2(it — |tq|)/efe and Q{x) 
the step function. The /c-dependent constant ZA,k is de- 
fined below {Za,a = !)• Since Rk=o{^) = 0, the action 
Sk=o coincides with the action (27) of the original model. 



For small /c, the function i^/c(q) gives a "mass" 



to 



the low-energy modes |q| ^ A: and acts as an infrared 
regulator as in the standard NPRG scheme [^(i, 17]. 
The scale-dependent effective action 

Tk[r.^] = -lnZfe[J*,J]+ / drV(j;(/>r + c.c.) 

-ASkir.^] (30) 

is defined as a (slightly modified) Legendre transform 
which includes the explicit subtraction of AS'/c[0*,^]. 
Here Z/e[J*, J] is the partition function obtained from 
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the action S + AS'/., Jj. a complex external source which 
couples linearly to the bosonic field V^r, and 



Mr) = (Mr)) 



S\nZk[J*,J] 

dJ*{T) 



(31) 



the superfiuid order parameter. The variation of the ef- 
fective action with k is governed by Wetterich's equa- 
tion [19], 

dkrkir.^] = ^TrjSfci^fc (ri'^[r,^]+^fc)"'}, (32) 

(2) 

where ^ is the second-order functional derivative of Ffc. 
In Fourier space, the trace in (32) involves a sum over 
momenta and frequencies as well as the two components 
of the complex field (/). 

We are primarily interested in two quantities. The first 
one is the effective potential defined by 



Vk{n) = —Tk[r,<l>] 



(33) 



where (/) is a constant (uniform and time-independent) 
field. The U(l) symmetry of the action implies that Vk{n) 
is a function of n = Its minimum determines the 

condensate density no,/c and the thermodynamic poten- 
tial (per site) Vo^k = yfc(^o,/c) in the equilibrium state. 
The second quantity of interest is the two-point vertex 



.(2) 
k,ij 



(r - r',r ■ 



(34) 



which determines the one-particle propagator Gk 



.(2)-l 



and enters the fiow equation (32). Here the in- 



dices i, j refer to the real and imaginary parts of < 
1 



Mr) 



V2 



[(/>ir(r) +^(/)2r(r)] . 



(35) 



Because of the U(l) symmetry of the action (27), the 
two-point vertex in a constant field takes the form [9] 



^(2) 
k,ij 



{q;(l)) = SijTA,k{q;n)^(i)i(i)jrB,k{q;n)^eijrc,k{q;n) 

(36) 

in Fourier space, where q — (q, zcj), is a Matsubara 
frequency and e^j the antisymmetric tensor. For q = 

(2) 

0, we can relate F^ to the derivative of the effective 
potential, 

rg,('/ = 0;« 

so that 



(37) 



rA,fe((/ = 0;n) = 14'(n), 

rBAq = 0;n) = vl'{n), 

rc,fe(g = 0;n) =0. 



(38) 



B. Initial conditions 

Since the action S + AS a = Sioc corresponds to the 
local limit, the initial value of the effective action reads 

Fa [r , (/>] = Floe [r ,</>]+/ drj^r (q)^q</>(q) , (39) 

^0 q 

where Fioc , ^] is the effective action in the local limit 
{t = 0). It is not possible to compute the functional 
Tioc 7 foi" arbitrary time-dependent fields [ ] . One 
can however easily obtain the effective potential V\oc{n) 

(2) 

and the two-point vertex F[^^ in a time-independent field 
(f). These quantities are sufficient to specify the initial 
conditions of the fiow within the approximations dis- 
cussed below. 

The initial effective action Fa [Eq. (39)] treats the local 
fiuctuations exactly but includes the intersite hopping 
term at the mean-field level, thus reproducing the strong- 
coupling random-phase approximation (RPA) [20-24]. 



C. Approximate solutions of the NPRG equations 

To solve the NPRG equations we expand the effective 
potential about its minimum. 



where 



Vo^k + 4(^- 



^0,/c) + yi^ - 



4 = 



dn 



At = 



Oil? 



no,fe)2, (40) 



(41) 



Vb,/c determines the thermodynamic potential per site 
and in turn the pressure P{fi,T) = — Vb,/c=o- We also 
use a derivative expansion where the two-point vertex 



r^^\q) = r)^\q;no k) in the equilibrium state is defined 

by 

^aMq) = ZaA^ci + 2c?i) + ^A.fcw^ + 4, 
rB,fe(q) = Afe, (42) 
rc,fe('2) = Zc,k(^- 
The initial values Za,k, Va,a, Zc,k, (5a and Aa are de- 
duced from F^^ and VA(n). The fiow equations of Z^,/c, 
VA.k^ Zc,k^ ^k and \k are then obtained from the exact 
fiow equation (32) within a simplified Blaizot-Mendez- 
Galain-Wschebor scheme [^5, 26]. We refer to Ref. [o] 
for a detailed discussion. 



.(2) 



D. Infrared behavior 

The infrared behavior can be obtained from the action 
TkW,4>] = dr j d^r[cp*{Zc,kdr - Vam^I 

- Z^.fetV' + 4)(i6 + Y(n - no.fc)' + Fo.fc] . (43) 
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Since we are interested in the low-energy limit, we con- 
sider the continuum limit where r becomes a continuous 
variable. Higher-order (in derivative or field) terms ne- 
glected in (43) give subleading contributions to the in- 
frared behavior. Most physical quantities of interest can 
be directly deduced from Eq. (43) [5]. The pressure is 
given by 



0.5 



(44) 



In the superfiuid phase, the superfiuid stiffness can be 
expressed as 



and the sound velocity reads 



c(/x,T) = 



k{ij,,T) ' 



(45) 



(46) 



where Hi = d^P/dfi^ is obtained from (44). 



IV. UNIVERSAL THERMODYNAMICS NEAR 
THE MOTT TRANSITION 

In this section, we first discuss the QCP between the 
superfiuid phase and the Mott insulator. We show that 
elementary excitations are quasi-particles with spectral 
weight Zqp, effective mass m* and effective "scattering 
length" a*. Zqp, m* and a* are computed as a function of 
t/U using the NPRG equations. We then verify that near 
the generic Mott transition thermodynamics quantities, 
as obtained from the NPRG approach, can be expressed 
in terms of the universal scaling functions introduced in 
Sec. II. 



A. Quantum critical point 

At the quantum critical point {fi = fic^T = 0) between 
the superfiuid and Mott insulating phases, the effective 
action F = Tk=o [Eq. (43)] takes the form 



r[0*,0] 



b*{Zcdr-ZAtV^ 



(47) 

up to a constant term f3NVo and neglecting higher-order 
(in field and derivative) terms. Equation (47) is valid at a 
generic QCP where the transition is driven by a density 
change. At a multicritical point, where the transition 
takes place at fixed (commensurate) density, Zq vanishes 
and the leading time-derivative term — Va^t^ must be 
included in the effective action [4, 5]. 

From Eq. (47), we can identify the elementary excita- 
tions at the QCP. On the lower part of the transition line 
(for a given Mott lobe), Zc is negative and it is conve- 
nient to perform a particle-hole transformation <p ^ (j)^ 



-0.5 



+ 




y 











/I 
/ I 



1.04 



1.02 



0.015 0.02 0.025 0.03 0.035 0.015 0.02 0.025 0.03 0.035 

ijU t/U 

FIG. 2. (Color online) Zc, A and Za vs t/U at the QCP 
between the superfiuid phase and the Mott insulator ft — 1. 
The + and — signs refer to the upper and lower parts of the 
transition line. 



(which changes the sign of the dr term in (47)). We can 
then define a quasi-particle field 



(r,T) = v1^</'(r,T), 



and rewrite the effective action as 



where 



dr 



2m* 



147ra* 
2 m* 



m 



m- 



\Zc\ 



\Zc\ 



(48) 



(49) 



(50) 



with m = l/2t the effective mass of the free bosons mov- 
ing in the lattice. We deduce from Eqs. (48,49) that 
elementary excitations are quasi-particles with mass m* 
and spectral weight 



\zcr. 



(51) 



They are particle-like if Zc > and hole-like if Zc < 0. 
The effective interaction between two quasi-particles is 
determined by the effective "scattering length" a*. 

The quantum phase transition at /i = —2dt between 
the superfiuid phase and the vacuum (which can be seen 
as a Mott insulator with vanishing density) is a partic- 
ular case of a superfiuid-Mott-insulator transition which 
differs from the superfiuid- vacuum transition discussed in 
Sec. II A only by the presence of the lattice. In this case, 
Za = Zc = I (the single-particle propagator is not renor- 
malized [ , ]), so that Zqp = 1 and m* = m = l/2t. Fur- 
thermore, the interaction constant A = 47ra/m = Sirta 
can be calculated analytically and related to the scatter- 
ing length 



1 



S7r{t/U + A) ' 



; 0.1264 



(52) 



of the bosons moving in the lattice [ ], which gives a* = a. 

For a generic QCP between the superfiuid phase and 
a Mott insulating phase with nonzero density (n = 
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FIG. 3. (Color online) Quasi-particle weight Zqp, effective 
mass m* and scattering length a* vst/tc at the QCP between 
the superfluid phase and the Mott insulator n = 1 (tc is the 
value of t at the tip of the Mott lobe). The QMC data are 
taken from Ref. [27]. In the bottom figure, the scattering 
length a of the free bosons in the lattice is given by Eq. (52). 
The + and — signs refer to the upper and lower parts of the 
transition line. 



1, 2, 3, • • • ), the values of Zq 



and a* can be obtained 



from the numerical solution of the NPRG equations. Fig- 
ures 2 and 3 show Za, Zc, A and Zqp, m*, a* as a function 
of t/U for the transition between the superfluid phase 
and the Mott insulator with density n = 1. The van- 
ishing of Zc at the multicritical point implies that m* 
vanishes while Zqp and a* diverge as we approach the tip 
of the Mott lobe. In addition to the NPRG results, in 
Fig. 3 we show m* obtained from quantum Monte Carlo 
(QMC) simulations [ ] as well as m* and Zqp obtained 
from the strong-coupling random-phase approximation 
(RPA) (see Appendix B). The RPA value for the hop- 
ping amplitude at the tip of the Mott lobe, tc/U o:^ 0.286, 
is far away from the QMC {tc/U = 0.034083) [ ] or 
NPRG {tc/U 0.0339) results [28]. Nevertheless the 
RPA predictions for the quasi-particle weight Zqp and 



FIG. 4. (Color online) Same as Fig. 3 but for the transition 
from the superfluid phase to the Mott insulator with density 
n = 2. 



the effective mass m*, when plotted as a function of t/tc, 
are in good agreement with the NPRG and QMC re- 
sults (Fig. 3). As expected, the results for the lower 
part of the transition line are trivial in the limit t ^ 0: 
Zqp = m* /m — a* /a = 1 (they are simply obtained by 
considering the motion of a hole in a Mott insulator with 
one boson per site). 

In Fig. 4 we show Zqp, m* and a* as a function of 
t/tc for the transition between the superfluid phase and 
the Mott insulator with density n = 2. The results are 
similar to the case of the transition to the first Mott lobe 
(n = 1) but the behavior for t ^ is different. The limit- 
ing values of the effective mass and quasi-particle weight 
are given by the strong-coupling RPA (see Appendix B), 



lim Zq 

t/c/^o 



lim 



m ( 

m* \ 



1 



(lower branch), 
(upper branch). 



(53) 

where n denotes the (commensurate) density of the Mott 
insulator and we distinguish between the upper {ji c^Un) 
and lower {ji U{n — 1)) branches of the transition line. 
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Since rrf /m and a* /a are typically of order one (except 
close to the Mott lobe tip), the characteristic energy scale 
l/m*a*^ below which the physics is universal is roughly 
set by the hopping amplitude t. As we approach the 
tip of the Mott lobe (point C in Fig. 1), m*a*'^ diverges 
and the energy scale l/m*a*^ vanishes. The low-energy 
physics is then controlled by the multicritical point. 



B. Universal thermodynamics 

Since the superfluid-Mott-insulator transition belongs 
to the dilute Bose gas universality class, near the QCP 
the thermodynamics can be expressed in terms of the uni- 
versal scaling functions introduced in Sec. II as well as the 
nonuniversal parameters m* and a* . To ensure that there 
is no other nonuniversal parameter, we must verify that 
the chemical potential (or, more precisely, dji = fi — jHc) 
couples to the elementary excitations with no additional 
renormalization. Slightly away from the QCP, the shift 
^/i = /i — /ic in chemical potential implies a change 



where 



SS 



r/3 

/ dr Y,rr^r (54) 



in the action. To lowest order in ^/i, SS induces a cor- 
rection 

ST[^\(t)] = -Z^6iJi dr j d^r(t)''(t) (55) 

to the effective action (47) at the QCP, where is a 
renormalization factor. Using the Ward identity = 
Zc (see Appendix C), we obtain 



= -sgn(Zc)^/i^ dr j d' 



r(l)*(l). (56) 



We conclude that sgn(Zc)(^/i acts as a chemical potential 
for the elementary excitations at the QCP [ ]. This 
implies that ±(5/i/T will enter scaling functions with no 
additional scale factor. This result agrees with general 
considerations on the scaling of conserved densities near 
a continuous quantum phase transition [ ]. 

Near the superfluid-Mott-insulator transition we can 
therefore write the pressure as 



(57) 



or 



3/2 



^(e) 87rV'2m*a*2|e| 



(59) 



is obtained from Eq. (8) by replacing m and a by m* 
and a*. The +/— sign in Eqs. (57,58) corresponds to 
particle/hole doping of the Mott insulator {i.e. the up- 
per/lower part of the transition line). Universality argu- 
ments imply that the singular part of the pressure can 
be expressed in terms of the scaling function but do 
not allows us to determine the regular part. To obtain 
the latter, we note that the compressibility hz = d'^P/dji^ 
vanishes in the T = Mott insulator and has therefore 
no regular part. 



(60) 



(see Eq. (9)). Integrating this equation twice with respect 
to /i yields Eq. (57) with Pc the value of the pressure at 
the QCP and fic the density at the QCP. 

The results obtained from a numerical solution of the 
NPRG equations show that near the superfluid-Mott- 
insulator transition the pressure can be expressed in 
terms of the universal scaling function T discussed in 
Sec. II. In the dilute classical regime, sgn(Zc)J/i < and 
^ T, our results are compatible with the expected 
result 



2^ 



3/2 



2^5/2g-|5/i|/T 



(61) 



(see Eq. (15)). However, a precise comparison is pre- 
vented by numerical difficulties due to the extremely 
small values of the pressure in this regime. The finite- 
temperature pressure at point B in Fig. 1 is shown in 
Fig. 5. We find a very good agreement with 



P(/Xc,T) = Pe + C(5/2) 



m 
2^ 



3/2 



2^5/2 ^g2) 



below a crossover temperature scale T ^ 2t. Figure 6 
shows the T = pressure at fixed t/U and for a density 
n varying between and 1 (see the dotted line in Fig. 1). 
Near the Mott insulating phases (n or n 1), for 
\5jii\ < t, there is a very good agreement between the 
NPRG result and the universal form 



P{fl,0) = Pc^UcSfl 



64 



1- — ^Vm*a*2|(5/i| . (63) 
157r 



(58) 



It should be noted that the agreement is better with the 
"Lee-Huang-Yang correction" (last term of (63)) than 
without. Differentiating (63) with respect to /i, we obtain 

n(M,0) = ne + ^ (^l-^x/m*a*2|5^|^, (64) 
- — \/m*a*^|J/i| ) . (65) 



^(m,o) 
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FIG. 5. (Color online) Pressure P(^\ic^T) vs temperature T 
[t/[/ = 0.02 and [icjV ^ 0.15 = (point B in Fig. 1)]. The 
dashed (green) line corresponds to Eq. (62). The inset shows 
a log-log plot and the T^^^ dependence of P(/ic, T) — Pc. 
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FIG. 6. (Color online) Zero-temperature pressure P(/i, 0) vs 
chemical potential /i along the dotted line in Fig. 1 (the inset 
shows the density n = dP/d/n). The bottom figures show the 
behavior near the Mott insulating phases n = and n = 1. 
The dashed (green) line corresponds to Eq. (63) and the dash- 
dotted (blue) one to the "mean-field" result P = Pc + nc^/i + 
((5//^)m*/87ra*. 
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FIG. 7. (Color online) Condensate density no(/i, 0) vs /i along 
the dotted line in Fig. 1. The insets show the behavior near 
the Mott insulating phases n = and n = 1. The dashed 
(green) line corresponds to Eq. (67) and the dash-dotted 
(blue) one to the "mean-field" result no = Zqpm* |^/i|/47ra*. 
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FIG. 8. (Color online) Superfluid stiffness ps(/i, 0) vs /i along 
the dotted line in Fig. 1. The insets show the behavior near 
the Mott insulating phases n = and n = 1. The dashed 
(green) line corresponds to Eq. (69) and the dash-dotted 
(blue) one to the "mean- field" result |^/i|/47ra*. 



perature, 



The condensate density no(/i,T) in the superfluid 
phase can be expressed in terms of the scaling function X 
[Eq. (12)]. However, since only the coherent part of the 
excitations {i.e. the quasi-particles) condenses, Eq. (12) 
can be used for the condensate density of the quasi- 
particles while I^P = Zqpl^p, which leads to 

no(M,r) = (^'I}^y\(^^^,g{5^^)^ (66) 

near the superfluid-Mott-insulator transition. The fact 
that no(/i,T), contrary to other physical quantities dis- 
cussed in this section, depends on the quasi-particle 
weight can be understood by noting that (/) is not invari- 
ant in the local gauge transformation (C3) and therefore 
not "protected" by the Ward identity (C4). At zero tem- 



no(M,0) = ^qp^^B^ " ^Vm*a*^M^ , (67) 

where m*|(5/i|/47ra* \n — nd is the density of excess 
particles (or holes) with respect to the commensurate 
density fic of the Mott insulator. The T = condensate 
density along the dotted line in Fig. 1 is shown in Fig. 7. 
Near the Mott insulating phases nc = and nc = 1, we 
find a very good agreement with Eq. (67). 

The superfluid stiffness can be expressed using the scal- 
ing function J [Eq. (13)], 



Psi,,T) = {^^) (68) 
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FIG. 9. (Color online) Sound mode velocity c(/i, 0) along 
the dotted line in Fig. 1. The dashed (green) line shows the 
expected result near the Mott transition [Eq. (71)]. 



Using the results of Sec. (II B), we obtain 



i(m,0) 



^ (69) 

Hi 

again in very good agreement with the NPRG approach 
(Fig. 8). ps, contrary to no, is independent of the quasi- 
particle weight Zqp. This follows from the fact that 
the superfluid stiffness can be related to the current- 
current density correlation function [ ], which is a 
gauge-invariant quantity. As a result, the ratio between 
condensate density and superfluid stiffness. 



^o(m,Q) 
m*ps(M,0) 



(70) 



(to leading order in m*a*^|(5/i|), explicitely depends on 
the quasi-particle weight while it is equal to unity in a 
dilute Bose gas. 

From Eqs. (65) and (69), we deduce 



c(/i,0) 



A^(/i,0) 



(71) 



to leading order in m*a*^|^/i|. As the condensate den- 
sity no and the superfluid stiffness p^, c is related to the 
density of excess particles (or holes) |n — nd rather than 
the full density n. The NPRG results show that Eq. (71) 
is very well satisfied near the Mott transition (Fig. 9). 

The superfluid transition temperature is determined 
by the scaling function 



n{g{T,)), 



(72) 



where 1-L is the universal scaling function introduced in 
Sec. II A. Using (26), we then obtain 



m* V87rC(3/2)a* 



2/3 



27r / |n - nc\ \ 
m* V2C(3/2); 



2/3 



(73) 



(to leading order in m*a*^|^/i|) near the T = Mott 
transition, in very good agreement with the NPRG re- 
sults (Fig. 10). 
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FIG. 10. (Color online) Superfluid transition temperature Tc 
vs 11 along the dotted line in Fig. 1. The dashed (green) line 
shows the expected result near the T = Mott transition 
[Eq. (73)]. 



C. RG flows and approach to universality 

Figure 11 shows the flow of the coupling constants Zc^k 
and A/e in the zero-temperature superfluid phase near a 
QCP. Exactly at the QCP {dji = 0), one can clearly dis- 
tinguish two regimes: i) a high-energy (or short-distance) 
regime k ^ kx where lattice effects are important and the 
dimensionless coupling constant [5] 



Zc,kZA,kt 



(74) 



is large, ii) a weak-coupling ( "Bogoliubov" ) regime k < 
kx where A/c <C 1 and the flow is governed by the Gaus- 
sian fixed point A = 0: A/^, Zc^k and ZA^k are then 
nearly equal to their fixed-point values [Eqs. (50)] while 
Xk (X k vanishes in agreement with its scaling dimension 
[Xk] = 4: — d — z = —1 at the Gaussian fixed point {d = 3 
and z = 2). In the momentum regime |q| < kx^ the 
quasi-particles with mass m* and scattering length a* 
introduced in Sec. IV A are well defined and the physics 
becomes universal. The crossover scale kx between the 
two regimes is typically of the order of A = \^ {k~^ is 
equal to a few lattice spacings). 

Away from the QCP, chemical potential and tempera- 
ture introduce two new momentum scales, the "healing" 
scale 



kh = V'2m*|(^/i| 
and the thermal scale 



kr = V2m*T. 



(75) 



(76) 



Universality requires kh^kr <C kx. Since kx 
(except close to the tip of the Mott lobe) these conditions 
can be rewritten as 



m*a*^|(5/i| <C 1, 
Vm*a*2T < 1. 



(77) 
(78) 
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FIG. 11. RG flows of Afe (Aa ^ 37), Xk and Zc,k at the QCP 
(t/[/ = 0.02, /ic/t/ ^ 0.15, T = 0) (point B in Fig. 1) and in 
the nearby superfluid phase 6/j/U — — 10~^ (A^^ ai 
A and Zc stand for Afc=o and Zc,fc=o, respectively. 



'SF \ 



In the low-energy limit the system behaves as a gas 
of weakly-interacting quasi-particles if the dimensionless 
coupling constant 



A. 



(79) 



is small. The last result in (79) is obtained using kh <C 
which allows us to approximate Zc^Uh^ ^A,kh 
by their k = values. Since kx ~ ^ 1, uni- 
versality {kh <^ kx) implies weak coupling (/c/^a* <C 1). 
Using Eq. (64), the weak-coupling/universality condition 
kha* <C 1 [Eq. (77)] can be rewritten as 



v4 



n — He o 



*3 



< 1. 



(80) 



Equation (80) is similar to the usual condition for a bo- 
son gas to be dilute except that it involves the excess 
density of particles (or holes) \n — nd (with respect to 
the commensurate density of the Mott insulator) rather 
than the full density n of the fluid. 

For k ^ kh^ Xk and Zc^k depart from their fixed-point 
values at ^/i = (Fig. 11) and vanish logarithmically be- 
low a "Ginzburg" momentum scale kc which is exponen- 
tially small at weak coupling {Xk^ <C 1). In a dilute Bose 
gas, the Ginzburg scale manifests itself by the appear- 
ance of infrared divergences in the perturbation theory 
about the Bogoliubov approximation. Although these di- 
vergences cancel out for thermodynamic quantities, they 
do have a physical origin: they result from the coupling 
between longitudinal and transverse (phase) fluctuations 
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FIG. 12. Gharacteristic momentum scales in a dilute super- 
fluid, at the generic Mott transition (/i = /ic), and in the 
superfluid phase near the Mott transition. (In the dilute su- 
perfluid, is of the order of the scattering length a, while 
it is of the order of the inverse lattice spacing in the Bose- 
Hubbard model.) 



and lead to a divergence of the longitudinal susceptibility 
- a general phenomenon in systems with a broken con- 
tinuous symmetry [ ]. For k ^ kc, the RG flow crosses 
over to a "Goldstone" regime where the physics is domi- 
nated by phase fluctuations. We refer to Refs. [ , 33, 34] 
for a detailed discussion of the infrared behavior in the 
superfluid phase. The various regimes of the RG flow are 
summarized in Fig. 12. 



V. CONCLUSION 

We have presented a detailed study of the thermody- 
namics of a Bose gas near the generic (density-driven) 
Mott transition in the framework of the Bose-Hubbard 
model. In the critical regime, the physics is governed 
by weakly interacting quasi-particles with quasi-particle 
weight Zqp, effective mass m* and "scattering length" a*. 
Thermodynamic quantities can be expressed using the 
universal scaling functions of the dilute Bose gas univer- 
sality class. They are independent of the quasi-particle 
weight and the only nonuniversal parameters entering the 
scaling functions are m* and a*. A notable exception is 
the condensate density no, which is proportional to Zqp, 
thus allowing us to determine the quasi-particle weight 
from a thermodynamic measurement once m* and a* are 
known. The NPRG enables to compute Zqp, m* and a* 
as a function of t/U. We find that the strong-coupling 
RPA, although rather inaccurate to determine the phase 
diagram, gives rehable estimates of Zqp, m* and a* as a 
function of t/tc (with tc the value of the hopping ampli- 
tude at the tip of the Mott lobe). 
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The thermodynamics of a two-dimensional Bose gas in 
an optical lattice has recently been measured near the 
superfluid- vacuum transition [ ]. A similar experiment 
in a three-dimensional gas would allow us to test the 
universality class of the generic three-dimensional Mott 
transition and the predictions of the NPRG approach re- 
garding the values of Zqp, m* and a*. A measurement of 
the temperature dependence of the pressure in the quan- 
tum critical regime at = [Eq. (62)] would directly 
provide us with the value of the effective mass m* . Quite 
interestingly, m* strongly varies with both the ratio t/U 
and the commensurate value of the density in the Mott 
insulator [Eq. (53)]. Measuring the scattering length a* 
and the quasi-particle weight Zqp is more challenging as 
it would require to reach temperatures much smaller than 
the crossover temperature T ^ 2t below which the ther- 
modynamics becomes universal, which is not possible yet 
in actual experiments. 
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Appendix A: Perturbative calculation of scaling 
functions 



In this Appendix, we briefly review the perturbative 
calculation of the universal scaling functions of the three- 
dimensional dilute Bose gas universality class (Sec. II B). 

At low temperatures and positive chemical potential, 
the scaling functions can be obtained from a one-loop 
calculation (Bogoliubov theory). To organize the loop 
expansion [ ], we introduce a parameter / (which will 
eventually be set to 1) and consider the partition function 



(Al) 

in the presence of a complex external source J. The 
superfluid order parameter is deflned by 



^, , l S\nZ[J*,J] 15lnZ[J*,J] 



(A2) 



3r(J*.^ + CO.), 



I dJ*{r,T) ' I dJ{r,T) 

We now introduce the effective action 

r[cl)*,cl)] = -jlnZ[J*,J]+ dT J d 

(A3) 

defined as the Legendre transform of the thermodynamic 
potential lnZ[J*, J]. The loop expansion is an ex- 
pansion in l/l. To one- loop order. 



r[0*,0] = s[r,4>] + ^Tvin(-g-i[0*,</']) + o(r2), 

(A4) 



where Qc [^* , ^] is the classical propagator 



^c[^,^';0*,0] = 



(A5) 



(we use the notation x = (r, r) and = (r',r')). 

All thermodynamic quantities can be obtained from 
the effective potential defined by 



(A6) 



{V denotes the volume of the system), where is a con- 
stant (uniform and time-independent) field. The U(l) 
symmetry of the action (27) implies that V{n) is a func- 
tion of the condensate density n = The minimum 
of the effective potential determines the condensate den- 
sity no and the thermodynamic potential Vq = V{no) 
per unit volume in the equilibrium state. The pressure 
is then given by 



P{^,T) = -Vo. 



(A7) 



To compute the effective potential to one-loop order, we 
need to evaluate the trace in Eq. (A4) with the classical 
propagator Qc evaluated in a constant field. 



-lUJri 



C-q-25|'^|' 

(A8) 

where ujn = n27r/P (n integer) is a bosonic Matsubara 
frequency, = Cq — /i and Cq = q^/2m. 



1. Zero temperature 

Let us first consider the zero-temperature limit. Per- 
forming the trace over Matsubara frequencies [ ], we 
obtain 



V(n) = -/^^+f + ^ 



/ (^q-^q-2^n)+(9(r2), (A9) 

Jci 



where f = f d^q/{27r)^ and 



(AlO) 



The condensate density no in the equilibrium state is 
obtained from F'(no) = 0, i.e. 



no 



9 21 



1 

q L^q 



(2^q + 3^no)-2 



•o(r^), (All) 



where £^q is defined by (AlO) with n = uq. 

Setting no = fi/g in the 0{l/l) term, we obtain the 
effective potential 



Ii)+0{1-'). (A12) 
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where = \/eq(eq + 2/i). To eliminate the dependence 
on the momentum cutoff A, we introduce the s-wave scat- 
tering length a, 



m 
Ana 

and rewrite Vq as 



9 Ja 



e(A-|q|) 

2ea 



(A13) 



Sna 
1 
2 



4 



(A14) 



setting / = 1. The sums over q are now convergent and 
we can take the infinite cutoff limit A ^ oo with a fixed. 
Using 



1 



1 



we finally obtained 



2 










^ ^3/2^5/2 


2i?J 


' ~ 157r2 



(A15) 



Vo 



64 
15^ 



\/ mo? II 



(A16) 



and in turn Eq. (17). Equations (18,19) are then deduced 
from n = —dVo/djii and /t: = dn/dja. 
Similarly, from (A11,A13) we deduce 



no 



4^~ 2 



to one- loop order (setting / = 1), which leads to Eq. (20) 
using 



1 



37r2 



(m/i)^/^ 



(A18) 



/q \ ^q 

The superfluid density Ug is defined by the variation 

ST = J d'riyef (A19) 

of the effective action when the superfluid order param- 
eter 0(r) = y/riQe^^^^^ acquires a phase slowly-varying in 
space. For a dilute Bose gas, Galilean invariance im- 
plies that Ug is equal to the fluid density, i.e. ns{/J^^O) = 
n(/i, 0). To leading order, the sound mode velocity 
c = psi K = ^Jns/m^^ is equal to ^J (i/m^ in agreement 
with the small-q behavior of E^^ = Y^eq(eq + /i). 

2. Finite temperature 

Similarly, we can compute the pressure P(0, T) at van- 
ishing chemical potential. For /i = 0, the condensate 



density no = in the equilibrium state. The effective 
potential Vq to one-loop order is then simply given by 
the non-interacting result. 



= -C(5/2) ( 
where ((x) is the Riemann zeta function 



(A20) 



3. Transition line 

The one-loop approximation fails near the superfluid 
transition temperature Tc [ ]. The transition tempera- 
ture can nevertheless be determined from a perturbative 
approach in the normal state, by considering the self- 
consistent one-loop self-energy correction (self-consistent 
Hartree-Fock approximation). 



2g 



E/- 



^9 I riB{e^ ^T.- ii), 

J q 



(A21) 



where nsix) = (e^^ — 1)""*^ is the Bose-Einstein distribu- 
tion function. The transition occurs when the renormal- 
ized chemical /i — S vanishes. 



f^ = 2g [ nB(eq) =25C(3/2) 

-J q 



27r 



3/2 



(A22) 



To lowest order g = 47ra/m, which gives Eq. (25). 

Appendix B: Strong-coupling RPA 

In the strong-coupling RPA [20-24], the effective ac- 
tion is given by (39). This implies that in the Mott in- 
sulating phase, the single-particle propagator takes the 
form 



where 



G(q, iuj) 



1 



Gioc(icj) 

-^qGloc(^^)' 



iuj -\- p — Un iuj -\- p — U{n — 1) 



(Bl) 



(B2) 



is the local propagator, n denotes the mean (integer) 
number of bosons per site. The instability of the Mott 
insulator is signaled by the appearance of a pole in the 
propagator, i.e. 1 — tq^oG^ioc(^^ = 0) = 0? which repro- 
duces the mean- field phase diagram [ ] . The value of the 
hopping amplitude at the tip of the Mott lobe is given 
by 



2dtc 



2n + 1 - 2Vn2" 



n. 



(B3) 
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On the transition line, we obtain 



reads 



G(q, iuj) 



for q, ^ 0, with 



m 



iuj — q^/2m* 



Gioc(O) 



(B4) 



(B5) 



where G[^^{iuj) = di^j G \oc {ioj) . As in Sec. IV A, we have 
performed a particle-hole transformation when we con- 
sider the lower branch of the transition line (hence the 
absolute value in Eq. (B5)). 

In the limit t ^ 0, Eq. (B5) leads to (53), a result 
which can be understood as follows. Let us add a particle 
at site r to a Mott insulator with n particles per site. In 
the limit t ^ 0, the only possible dynamics is due to the 
motion of the additional particle. Hopping of this particle 
between sites r and involves the matrix element 

(n;r|(g)(n + l;r'|t?/iJ,?/;r|n + l;r)(g)|n;rO =t(n + l) (B6) 

if we denote by (8)rJ^i;i*i) the state with rii particles at 
site r^. The "particle" eigenstates are therefore plane 
wave states. 



1^) = ^E^'^"l^ + 1'^) ^rvr (B7) 



with a dispersion law Aq = Aq — 2t(n + 1) Xli=i ^osg^ 
(with Aq a constant which takes the value 2dt(n + 1) at 
the quantum critical point), which leads to an effective 
mass vff jm = l/(n -|- 1). The single-particle propagator 



G(q,icj) = 



l(q|^nq)|Q)l^ 

iUJ — Aq 



(B8) 



where |0) = (8)r|^; r) denotes the ground state of the Mott 
insulator without the additional particle (in the limit t 
0). We deduce the quasi-particle weight 



Zqp = Kq|^/;t(q)|0)|2 = n+l. 



(B9) 



A similar reasoning for the motion of a hole leads to 
^qp = m/m* = n. 



Appendix C: Ward identity — Zc 

To lowest order in the effective action (47) at the 
QCP is modified by 



5T[(t)\(t)] = -Z^Sii dr j d^r(t)''(t). 



(CI) 



This implies that the effective potential is given by 

V{^,n) = V{^c.n)-Z^n5^. (C2) 
The invariance of the action S in the local (time- 
dependent) gauge transformation 

V^r ^ V^re'"", V^r ^ V^Je"^", /i ^ /i + zS^a (C3) 

implies that Zc = Zcif^c) satisfies the Ward identity [5] 



Zc 



djidn 



IJ,c,n=0 



which gives 



Zc = Zn. 



(C4) 



(C5) 
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